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The local composite operator A 2 ^ is added to the Zwanziger action, which implements the 
restriction to the Gribov region (l in Euclidean Yang-Mills theories in the Landau gauge. We 
prove that Zwanziger's action with the inclusion of the operator A 2 ^ is renormalizable to all 
orders of perturbation theory, obeying the renormalization group equations. This allows to 



o 



c3 



Oh, study the dimension two gluon condensate {A 2 ^ by the local composite operator formalism 

when the restriction to the Gribov region Q is taken into account. The resulting effective 
action is evaluated at one-loop order in the MS scheme. We obtain explicit values for the 
Gribov parameter and for the mass parameter due to {A 2 }, but the expansion parameter 
turns out to be rather large. Furthermore, an optimization of the perturbative expansion in 
order to reduce the dependence on the renormalization scheme is performed. The properties 
of the vacuum energy, with or without the inclusion of the condensate {A 2 }, are investigated. 
In particular, it is shown that in the original Gribov-Zwanziger formulation, i.e. without 
the inclusion of the operator A 2 the resulting vacuum energy is always positive at one- 
loop order, independently from the choice of the renormalization scheme and scale. In the 
presence of \A 2 ), we are unable to come to a definite conclusion at the order considered. In 
the MS scheme, we still find a positive vacuum energy, again with a relatively large expansion 
parameter, but there are renormalization schemes in which the vacuum energy is negative, 
albeit the dependence on the scheme itself appears to be strong. Concerning the behaviour of 
the gluon and ghost propagators, we recover the well known consequences of the restriction to 
the Gribov region, and this in the presence of {A 2 }, i.e. an infrared suppression of the gluon 
propagator and an enhancement of the ghost propagator. Such a behaviour is in qualitative 
agreement with the results obtained from the studies of the Schwinger-Dyson equations and 
from lattice simulations. 
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1 Introduction. 



The dimension two condensate (-A^) has received a great deal of attention in the last few years, 
see for example 011131101313111111011111112111311111131131111. This condensate was 
already introduced in :18 a in order to analyze the gluon propagator within the Operator Product 
Expansion (OPE), while in 05] the condensate (A|) was considered in the Coulomb gauge. A 
renormalizable effective potential for \AV) has been constructed and evaluated in analytic form 
up to two-loop order in the Landau gauge within the local composite operator (LCO) formalism 
in [3 Q2| ■ The output of these investigations is that a non- vanishing condensate is favoured as 
it lowers the vacuum energy. The renormalizability of the local composite operator formalism, 
see [201 f° r an introduction to the method, was proven to all orders of perturbation theory, in 
the case of \AV), in IT] using the algebraic renormalization technique [25- Besides the Lan- 
dau gauge, the method was extended to other gauges as, for instance, the Curci-Ferrari gauge 
[221123]) the linear covariant gauges |24[ I25| and, more recently, the maximal Abelian gauge |26j . 

As a consequence of the existence of a non- vanishing condensate (At), a dynamical mass param- 
eter for the gluons can be generated in the gauge fixed Lagrangian, see |H1 1121 [23] . We mention 
that a gluon mass has been proven to be rather useful in the phenomenological context, see e.g. 
[271 1281 155] . Moreover, mass parameters are commonly used in the fitting formulas for the data 
obtained in lattice simulations, where the gluon propagator has been studied to a great extent 
in the Landau gauge [33 EJJ [33 [33 [31 [33 [33 . 

The lattice results so far obtained have provided firm evidence of the suppression of the gluon 
propagator in the infrared region, in the Landau gauge. Next to the gluon propagator, also the 
ghost propagator has been investigated numerically on the lattice [31 [33 133 133 13H] , exhibiting 
an infrared enhancement. It is worth remarking that, in agreement with lattice results, this 
infrared behaviour of the gluon as well as of the ghost propagator has been obtained in the 
analysis of the Schwinger-Dyson equations, see [33 Ell El E3 E3 El S3 E3; as well as i n a 
study making use of the exact renormalizaton group technique |47j . 

The aim of the present work is to investigate further the condensation of the operator A^ 
in the Landau gauge using the local composite operator formalism. This will be done by taking 
into account the nonperturbative effects related to the existence of the Gribov ambiguities |48j . 
which are known to affect the Landau gauge fixing condition, d^A^ = 0. As a consequence of the 
existence of the Gribov copies, the domain of integration in the path integral has to be restricted 
in a suitable way. Gribov's orginal proposal was to restrict the domain of integration to the 
region f2 whose boundary d£l is the first Gribov horizon, where the first vanishing eigenvalue 
of the Faddeev-Popov operator, —d^ [d^5 ab + gf acb A c ^j, appears [IE]. Within the region Q the 
Faddeev-Popov operator is positive definite, i.e. —d^ [d^5 ab + gf^Afy > 0. One of the main 
results of Gribov's work |3H] was that the gluon, respectively ghost propagator, got suppressed, 
respectively enhanced, in the infrared due to the restriction to the region f2. 

In two previous papers [Ml 133] , we have already worked out the consequences of the restric- 
tion to the Gribov region f2 when the dynamical generation of a gluon mass parameter due to 
Afy takes place, also finding an infrared suppression of the gluon and an enhancement of the 
ghost propagator. In [33, we closely followed the setup of Gribov's paper J3HJ. In this work, 
we shall rely on the Zwanziger local formulation of the Gribov horizon. In a series of papers 



2 



|51| I52|. Zwanziger has been able to implement the restriction to the Gribov region $7 through 
the introduction of a nonlocal horizon function appearing in the Boltzmann weight defining the 
Euclidean Yang- Mills measure. More precisely, according to |51| I52|. the starting Yang- Mills 
measure in the Landau gauge is given by 



dfiy = DASQdpA*) det(M)e-( SYM+ ^ H ) , (1.1) 



where 



M ah = -dJdJ ab +gf acb A c ll ) , (1.2) 



and 



H = J d 4 xh(x) =g 2 J d 4 xf abc A b ^ {M- l ) ad f dec A e ^ , (1.4) 

is the so-called horizon function, which implements the restriction to the Gribov region. Notice 
that H is nonlocal. The parameter 7, known as the Gribov parameter, has the dimension of a 
mass and is not free, being determined by the horizon condition 

(h(x))=A(N 2 -l) , (1.5) 

where the expectation value (h(x)) has to be evaluated with the measure d[i~. To the first order, 
the horizon condition (jl.5j) reads, in d dimensions, 

_ N (d-1) 2 f d d k 1 

~ 4 9 J (^FTW ' ( } 

This equation coincides with the original gap equation derived by Gribov for the parameter 7 



Albeit nonlocal, the horizon function H can be localized through the introduction of a suit- 
able set of additional fields. As shown in the resulting local action turns out to 
be renormalizable to all orders of perturbation theory. Remarkably, we shall be able to prove 
that this feature is preserved when the local operator A 2 ^ is introduced in the Zwanziger action. 
Moreover, the resulting theory turns out to obey a homogeneous renormalization group equation. 
These important properties will allow us to study the condensation of the operator A^ within a 
local renormalizable framework when the restriction to the Gribov region $7 is implemented. 

It is worth remarking that the Gribov region is not free from gauge copies [541 1551 1561 I57j . 
i.e. Gribov copies still exist inside f2. To avoid the presence of these additional copies, a further 
restriction to a smaller region A, known as the fundamental modular region, should be imple- 
mented. At present, a clear understanding of the role played by these additional copies appears 
to be a very difficult task. Nevertheless, we should mention that, recently, it has been argued 
in [^Sj that the additional copies existing inside could have no influence on the expectation 
values, so that averages calculated over A or f2 might give the same value. 

The paper is organized as follows. In section 2, we give a short account of how the nonlo- 
cal horizon functional H can be localized by means of the introduction of additional fields. In 
section 3, we prove the renormalizability, to all orders of perturbation theory, of Zwanziger's 
action in the presence of the operator A^, introduced through the local composite operator for- 
malism. As the model has a rich symmetry structure, translated into several Ward identities, 
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it turns out that only three independent renormalization factors are necessary. The resulting 
quantum effective action obeys a homogeneous renormalization group equation, as explicitly 
verified at one-loop order. From this effective action, two coupled gap equations, associated to 
the condensate {A 2 ^ and to the Gribov parameter 7, are derived. Section 4 is devoted to the 
study of these gap equations at one-loop order in the MS renormalization scheme. It is worth 
mentioning that, under certain conditions, we find that it is possible that the condensate \AV) 
is positive when the horizon condition is imposed. We prove that in the MS scheme, and at 
one-loop order, the solution of the gap equations is necessarily one with (A 2 ^ > 0. We recall 
that without the restriction to the Gribov region $7, the value found for (A 2 } using the local 
composite operator formalism is negative, see 1121 125 j . Let us also mention here that in 
[El El El EH j a positive estimate for (A 2 ^ was obtained when using the OPE in combination with 
\A 2 }. These works were based on the observation of a certain discrepancy at relatively large 
momentum between the expected perturbative behaviour and the obtained lattice behaviour of 
e.g. the effective strong coupling constant and gluon propagator. This discrepancy could be 
accounted for by power corrections in \, due to a positive (^)qpe ghion condensate. The 
presence of such power corrections has also been discussed in |59| . We do not know if there 
is a direct connection between the condensate (A 2 ^ that we determine, and ^t4^ ope , as the 
latter is expected to contain only infrared contributions, according to an OPE treatment, while 
the gap equations fixing the gluon condensate and the Gribov parameter are evaluated using 
perturbation theory, implying that reliable results are only to be expected at a sufficiently large 
scale. 

Although the expansion parameter proves to be rather large, an attempt to obtain explicit 
values for the Gribov and gluon mass parameter is still presented. Also, we shall prove that in 
the original Gribov- Zwanziger model, the vacuum energy is always positive at one-loop order, 
irrespective of the choice of renormalization scheme and scale. We outline the importance of the 
sign of the vacuum energy, as it is related to the gauge invariant gluon condensate (F 2 U ), via 
the trace anomaly. From 

a _ P(d ) p 2 7 x 
~~ 2g 2 ^ ' 

the vacuum energy can be traced back to the value of the gluon condensate \F 2 ). In particular, 
for N = 3, from this anomaly one deduces 

/ J^FlS = -^E vac , (1.8) 



\4vr 2 "V " 11 

where the one-loop /3-function has been used. Hence, a positive vacuum energy implies a negative 



value for the condensate ( -§^F 2 \. This is in contradiction with what is found. In QCD, with 



quarks present, one can extract phenomenological values for \j^F 2 u j via the sum rules pU] . 

obtaining positive values for this condensate. It was discussed in [EJ how to obtain an estimate 
for it by means of lattice calculations. In the case of N = 3 Yang-Mills theory without quarks, 
it was found that 

(hz F li) = °- 14 ± °- 02GeV4 • (!- 9 ) 

Let us mention here that the Yang-Mills /^-function is negative up to the (known) four-loop or- 
der |H2l EH1 El] • Hence, E v & c and (-^iF 2 J\ will continue to have opposite sign at higher order. 
From this viewpoint, it seems to us that it would be an asset that the vacuum energy obtained 
from any kind of calculation is at least negative. 
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In section 5 we present an optimized expansion in order to reduce the dependence on the choice 
of renormalization scheme to a single parameter bo, related to the coupling constant renormal- 
ization. This is achieved by exchanging the mass parameters by their renormalization scale and 
scheme invariant counterparts and by re-expanding the series in the one-loop coupling constant. 
For 60 = 0, which corresponds to the MS scheme, we find a positive (AV), positive E VSLC and 
hence negative {F®/}. However, we find that a region of bo is existing in which the vacuum 
energy is negative, but unfortunately the dependence on 60 m this region happens to be very 
large. A higher order analysis seems to be required to reach more definite conclusions about the 
sign of (A*), E vac or (F%). 

For the benefit of the reader, we provide in section 6 an overview of some important conse- 
quences stemming from the presence of the Gribov and gluon mass parameters on the gluon and 
ghost propagators. We point out a particular renormalization property of the Zwanziger action 
in order to ensure the enhancement of the ghost propagator. Conclusions are written down in 
section 7, while the technical details of our work have been collected in the Appendices A and 
B. 



2 Local action from the restriction to the Gribov region. 

As explained in |51[ I52|. the nonlocal functional H can be localized by means of the introduction 
of a suitable set of additional ghost fields. More precisely, for the localized version of the measure 
d/i 7 we get, 

dfij = DADbDcDcDtpD^DojDu]e~ s , (2.1) 

where S is given by 1 

S = So-7 2 gj d*x {f abc A^ + f abc Affi) , (2.2) 

while 

So = S YM + J d^x [b a d^ + c a d^ (D M c)°) 

+ J d 4 x (Tf%d v (d v ^ c + gf abm A b u ^ c ) - uf<d v (d v u™ + gf abm A b u u™) 

-g(d v u?;)r b ™(D u c) b <p™) . (2.3) 

The fields (jp^f, (p^f) are a pair of complex conjugate bosonic fields. Each field has 4 (A 2 — l) 2 
components. Similarly, the fields (p^f , ca^) are anticommuting. The local action (|2.2j) is 
renormalizable by power counting. More precisely, it has been shown in |51| I52| I53| that the 
Green functions obtained with the action So with the insertion of the local composite operators 
jabcj^a^bc anc j jabcj^a-^bc are renorm alizable, the action So being indeed renormalizable by a 
multiplicative renormalization of the coupling constant g and of the fields |51l 1521 153j . We re- 
mark that the action Sq displays a global U (/) symmetry, / = 4 (A 2 — l) , with respect to the 
composite index i = (/x,c) = 1, ...,/, of the additional fields (<^ c , <^ c ,tJ° c ,u;° c ). Setting 

{<Pn ,<Pit > w m ) = VPiifii^ii^i) > \ 2A ) 

1 Our conventions are different from those originally used by Zwanziger. These can be obtained from ours by 
setting ip — > — ip and u) — * —uo. 
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we get 



So = S YM + J d A x {b a d^A1+lfd^{D^) a ) 

+ J 'tfx(l%d v (D v( p i ) a -Z3Zd u (D v LJ i ) a - 9 («)r frm (^c)Vr)) • (2-5) 

For the U(f) invariance we have 
UijSo = , 

u " = J* + - Vik) ■ (2 ' 6> 

The presence of the global U(f) invariance means that one can make use of the composite index 
i = (jjl,c). By means of the diagonal operator Qf = Ua, the i- valued fields turn out to possess 
an additional quantum number. As shown in |511 1521 I5"5| . the action Sq is left invariant by the 
following nilpotent BRST transformations, 







r , 




sc a 


= \gf ahc 


, b c c 


1 


sc a 


= b\ 


sb a 


= o, 


stf 


= 


SLJ? 


= 0, 


sul? 






= 0, 



(2.7) 
with 

sS = . (2.8) 

For further use, the quantum numbers of all fields entering the action Sq are displayed in Table 
1. It is worth noticing that, when f abc A a ^ip b ^ and f^A^Tp 1 ^ are treated as composite operators, 
they are introduced in the starting action So coupled to local external sources M" 1 , V™, namely 

- J d'x {Mf (D^iT + Vf (D^f) . (2.9) 

The horizon condition (|1.5j) is thus obtained from the quantum action by requiring that, at the 
end of the computation, the sources M!f , V™ attain the physical values, obtained by setting 

Mf v = V$ = 7 2 S a % u . (2.10) 



Indeed, expression (|2.9j) reduces precisely to that of eq. (|2.2j) when the sources , V™ attain 
their physical value. 





Aa 
A P 




c a 












dimension 


1 





2 


2 


1 


1 


1 


1 


ghostnumber 





1 


-1 











1 


-1 


Q /-charge 














1 


-1 


1 


-1 



Table 1: Quantum numbers of the fields. 
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M» 


N ai 




77 


r 


K a 


L a 


dimension 


2 


2 


2 


2 


2 


2 


3 


4 


ghostnumber 


-1 





1 





-1 





-1 


-2 


Q /-charge 


-1 


-1 


1 


1 















Table 2: Quantum numbers of the sources. 



3 Renormalization of the Zwanziger action in the presence of 
the composite operator A® A®. 

The purpose of this section is to show that the renormalizability of the local action So is preserved 
when, besides the operators f abc A^ip'jf and f abc A a J^, also the local composite operator A^A® 
is introduced. This is a remarkable feature of the Zwanziger action, allowing us to discuss the 
condensation of the operator A a A a when the restriction to the Gribov region 0, is implemented. 
To discuss the renormalizability of the model in the presence of A 2 , we start from the following 
complete action 

S = So + S s + S ext , (3-1) 

where S s is the term containing all needed local composite operators with their respective local 
sources, and is given by 

S s = s J d'x (-UJ? (D^iT ~ V* (D^f - U™V™ + \r,A%Al - \crr^ , (3.2) 

where the BRST operator acts as 

sUf = M? 

„yai _ Ajai 

[A H 

and 

ST] = T , ST = . (3-4) 

Therefore, for S s one gets 

Ss = J d 4 x (-M« (D^r - gU™f abc (D,c) b + U« {D.uif 

- M™V™ + U«N$ + \rA-A- + rjA^c" - ^r 2 ) ■ (3.5) 

As already noticed, the sources M™, V? 1 are needed to introduce the composite operators 
(D^ipi) 01 and (D fJ Jp i ) a . The sources U® 1 , N® 1 define the BRST variations of these operators, 
given by {D^ojif and (D p Zoi) a . The physical value of these sources is given by 

K = < = 0- (3.6) 

The local composite operator A a ^A a ^ and its BRST variation, A^d^c 0- , are then introduced by 
means of the local sources r, rj. We also notice that the complete action £ contains terms 
quadratic in the external sources, namely (MffV™ - U™N™) and (t 2 . These terms, allowed 
by power counting, are in fact needed for the multiplicative renormalizability of the model. As 



= 0, 
0, 



(3.3) 
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shown in the dimensionless LCO parameter £ of the quadratic term in the source r is needed 
to account for the divergences present in the correlation function ( s A 2 l (x)A 2 ,(y)^ for x — > y. It 
should be remarked that, unlike for the term quadratic in the external source r, we have not 
introduced a new free parameter for the quadratic term (MffV?* — UffNjf^j in expression (|3.5|) . 
As we shall see, this term goes through the renormalization without the need of introducing a 
new parameter for its renormalizability. This is a remarkable feature of the Zwanziger action 
which plays an important role when the ghost propagator in the presence of the Gribov horizon 
will be discussed, see section 6. 

Finally, the term S* ex t is the source term needed to define the nonlinear BRST transformations 
of the gauge and ghost fields, i.e. 



(3.7) 



The technical details concerning the algebraic renormalization procedure have been worked out 
in Appendix A. In summary, the Zwanziger action in the presence of the local operator A^A^ 
is multiplicative renormalizable. In turn, this ensures that the quantum effective action obeys 
the homogeneous renormalization group equations (RGE). This is an important feature of the 
model, which will be useful when we shall try to obtain estimates for both the Gribov and mass 
parameter. 

The effective action is defined upon setting the sources Uj*, Nj£, K£, L a and 77 equal to zero 
and implementing the condition ([2. 10)1 . Doing so, we get 

S = S + S, + j d±x(^-AlAl- Q -T 2 

S, = |d 4 x[- 7 V afec ^<-7V abc ^ c -4(iV 2 -l)7 4 ] • (3.8) 

The term -4 (N 2 - l) 7 4 originates from the quadratic term in the external sources, namely 
(— Mj^V™ + U^N™), in expression (|3.5|) . evaluated at the physical values given by ea. Q2.10l) . 

Following [3J 1121 1201 125j , we introduce a Hubbard-Stratonovich field a by means of the following 
unity 

1 = f[da]e~^^ dAx \-^ Aa ^~ CT \ , (3.9) 



to remove the term proportional to r 2 . The source r is henceforth linearly coupled to the field 
cr, as can be directly seen from the action, which now reads 



S = Sq + 5 7 + S a + J d x y T ~ J ' 



The following identification is easily derived 1121 1201 125) 

1 



A%A*) = --(*) , (3.11) 



from which it follows that a nonvanishing vacuum expectation value of the field a will result in 
a nonvanishing condensate {A a Afy. 
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., - 0. (3.14) 



The quantum action T is obtained through the definition 

e~ r = J [d&] e ~ So ~ s '>- s ° , (3.12) 
where $ is a shorthanded notation for all the relevant fields. 

The value for (a) is found through the minimization condition 

dT 

sf=°- (3 - 13) 

The horizon is implemented by the condition |51l I52j . 

8T 

Let us show this here. The following equivalence is readily found 

|L = & (gf abc A^) + (gf abc Affi) = -8 {N 2 - 1) 7 2 , (3.15) 

^From expressions (|l.lj) and (|2.2|) . it follows that 

-2 7 2 (h) = (gf abc A a ^) + (gf^Affi) . (3.16) 

The combination of ea. ()3.15|) with ea. (|3.16|) gives rise to the horizon condition ea. (|1.5[) . In order 
to conclude this, it is tacitly assumed that 7 / 0. We notice that the condition ()3.14j) does 
possess the solution 7 = 0. This is an artefact of the reformulation of the horizon condition in 
terms of the equation (|3.14j) . and must be excluded as it does not lead to the horizon condition 
(jl.5j) . We shall, however, continue to keep this solution of the gap equation (|3.14j) . as 7 = 
corresponds to the case where the restriction to the Gribov region £1 would not be implemented. 
In this case, we must only solve the gap equation stemming from ea. ()3.13|) with 7 = 0. 

The original Gribov- Zwanziger model, i.e. without the inclusion of the operator A 2 ^ is ob- 
tained by only retaining the condition (|3.14j) with a = 0. 

Up to now, the LCO parameter £ is still a free parameter of the theory. We do not intend 
here to give a complete overview of the LCO formalism, we suffice by saying that £ is fixed by 
the demand that the action T should obey the homogeneous renormalization group equation 

+ v {g2) w + lAg2)l2 w + >(5> ^) r = ' (3 - 17) 



with 



— w ! I 2\ 2 

= 772(5 h 



dp 

T^z = ia{g 2 )(T. (3.18) 

This can be accomodated for by making Q a function of the running coupling constant g 2 , in 
which case it is found that 

C(g 2 ) = ^ + <i + (29 2 + --- ■ (3.19) 

We refer to the available literature QH H3 |H H3 EH for a detailed account of the LCO 
formalism. 
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3.1 Renormalization group invariance of the one-loop effective action in the 

MS scheme without the inclusion of A 2 ,. 

a* 

Before proceeding with the detailed analysis of the horizon condition in the presence of the local 
operator A^A^, let us first derive the horizon condition and check the explicit renormalization 
group invariance of the quantum action T by switching off the source r coupled to the operator 
A^A^. This amounts to consider the original Gribov-Zwanziger model. We consider thus the 
action 

S = S + 5 7 . (3.20) 
The one-loop effective action is easily obtained from the quadratic part of ea, ()3.20|) 

-r (1) / [nr- -S' 



[£>$] e r b ^ d , (3.21) 



with given by 

Squad = I (1 X 



\ (d,A a u - d u Alf + i- {d.Alf + 



I 2 



g (f abc A^ + r bc Afi%) - 4(iV 2 - 1) 7 4 



(3.22) 



where the limit a — > is understood in order to recover the Landau gauge. After a straightfor- 
ward computation, one gets 

rW = -4(iV 2 - 1) 7 4 + (d - 1) | In (p 4 + 2A^V) . (3.23) 

Dimensional regularization, with d = 4 — e, will be employed throughout this work. Taking the 
derivative of r' 1 ', one reobtains the original gap equation for the Gribov parameter 7 , namely 

dr« N(d-1) o f d d p 1 

-0=^1 = — i -g 2 ^r— : 5—^. 3.24 



5 7 4 " J (27r) d (p 4 + 2iV 5 2 7 4 ) 

More precisely, recalling that 

ln(p +p ) =_ — In -3-3 + -—J, (3.25) 



( 27 r) d V ' ~ 32vr 2 V 7i 4 J e 32vr 2 ' 

the one-loop effective action reads 



r »> = - i(N * _ lh . _ ^ ( W ( ln ^ - |) , (3.26) 
where the MS renormalization scheme has been used. 

In order to check the renormalization group invariance of r« we need to know the anoma- 
lous dimension of the Gribov parameter 7. This is easily obtained from ea. (|A.34|l . yielding 

lA9 2 ) = -\(^-lA{9 2 )] , (3.27) 



' 7 vy ; 2 V 2 5 2 
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where 7^ (g 2 ) stands for the anomalous dimension of the gauge field A ^. Thus, at one-loop order, 
Furthermore, from (see e.g. |65j ^ 

it follows 

p— =0, (3.30) 
which establishes the RGE invariance of the effective action at the order considered. 

We are now ready to face the more complex case in which the local composite operator A^A^ 
is present. This will be the topic of the next section. 



4 One-loop effective action in the MS scheme with the inclusion 

of Al. 

4.1 Calculation of the one-loop effective potential. 

Let us turn to the explicit one- loop evaluation of the effective action T in the presence of ^4 2 . 
At one- loop, it turns out that 2 



-4(iV 2 -l) 7 4 + 



a 



2gK 



+ 



N 2 -\ 



In det 



' 2 2iV g V 

p <v 5 — 



'V'- - l>i>l>" ( 1 - - ) + 37^, 



(4 



1) 



or 



d d p 
(2^ 



In 



4 1 o AT 2 4 1 5 cr 2 



(4.2) 



Before calculating the integral, we quote the two gap equations 



dT o ,d-l ,„ 

- = ^ 7 3 = 7 3 — ^ 



(iV 2 - 1)^-1) f d 



2 



Co 



(2vr) V + f + 2iV 5 2 7 4 



0, 



(2tt) V + f P 2 + 



(4.3) 



The second gap equation of (|4.3j) . being the horizon condition, gives rise to the one obtained in 
the previous paper [IHj > while the first one describes the condensation of ^4 2 when the restriction 
to the Gribov region Q is implemented. We notice that that the explicit value of the Gribov 
parameter 7 is influenced by the presence of (^4 2 )- 



It remains to calculate 



/(3>[ ?4+i * 2+c1 ' 



(4.4) 



2 We shall drop from now on the superscript indicating that we are working at one-loop order. 
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with 

Since 
with 

one has 



9<? 



2„,4 



c = 2A/y 7 



Co ' 

p A + bp 2 + c = (p 2 + wi) (p 2 + W2) , 
b+Vb 2 - 4c b-Vb 2 - 4c 



W2 



(2VT) 



^ In (p 2 + u 2 ) 



(4.5) 
(4.6) 
(4.7) 
(4.8) 



To make sense, the expression (|4.4|) should be real to ensure that the one-loop effective action 
is real-valued. Therefore, we must demand that c > 0. If b > 0, X is certainly real. However, 
when b 2 — 4c < 0, then also b < is allowed. We should thus have a positive Gribov parameter 
7 4 , while the condensate (^4 2 ) can be negative or positive, depending on the case. 



Using 



it holds 



d p _ / n 2\ 
-7 in [p + to 

(27r) d V ' 



— TO 

32^2 



2 , m 2 3 
-" ln - +2 



Ll) 



2_ ui 3 
32tt 2 U ~ p 2 + 2 



32vT 2 



2 w 2 3 

- - In =2 + o 
e /r 2 



Finally, in the MS scheme, we obtain 



r 



-4(J\T 



1 )7 4 + ^(l"^ 



2Co 



Co" 



+ 



3 N 2 



( 22. 

u 



+ 



or 



1287T 2 



In 



9£_ 

Co 



+ 



8^ 2 A^7 4 



2Ti 2 



+ 



Co 



Co 2 



8g 2 NY 



128tH 



hi 



2£ 

Co 



C 2 



2/I 2 



To lighten the notation a bit, let us introduce the new variables 3 



TO 



8# 2 iV7 4 

Co ' 



in which case the action (|4.11j) can be rewritten as 

r = - v _ —: — + 



[N* - 1) A 4 Com 4 



+ 



2g 2 N 
3 (N 2 - 1) 
256?r 2 



2<? 2 

^TO 2 + 



A 4 



In- 



??•?■' 



+ v 7 ' 



77V 



A 4 



+ ^to 2 — \/ m 4 — A 4 ^ 



hi 



27I 2 

m 2 - Vm 4 - A 4 5\ 
~ 6 



■772 



27i 



(4.9) 
(4.10) 



(4.11) 



(4.12) 
(4.13) 



(4.14) 



In comparison with the previous article |49| . we have the correspondence A —4"/ with the Gribov parameter 
7 4 as defined there. It is actually this 7 4 which will enter the modified propagators, see |49| and further in this 
paper. 
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We notice that the foregoing expression is also valid, i.e. real-valued, in the case in which 
m 4 < A 4 , as £+(m, A) and £-.(m, A), defined by, 



£+(m,\) 

£_(m,A) 
are complex conjugate 4 . 



m 2 + \Jm 4 - A 4 



In 



m 2 + a/ 



A 4 



6 



A 4 



In- 



27Z 2 

m 2 - Vm 4 - A 4 5 
6 



27i 2 



(4.15) 



The horizon condition, ea. (|3.14|) . can be translated to 



ar 

dX=° 



and the gap equation (J3.13|) to 



or 

dm 2 







m 



(4.16) 



(4.17) 



reduces to the result 



As a check of this one-loop calculation, the expression (|4.14j) with 
obtained earlier in ea, (|3.26j) . i.e. the original Gribov-Zwanziger model without the inclusion of 
A 2 . If A = 0, i.e. the case where the condensation of A 2 is investigated without implementing 
the restriction to the Gribov region 0, ea. (|4.14|) coincides with the result of [31 1121 125] . From 
|llj . one knows that 



^(9 2 )(Al) 



2g 2 



+ lA{g 2 )){Al) 



or, using the relation (|3.11|) and the definition ()4.13j) . 



_dm 2 



Tm 2 ' 



2g 2 



7A{g 2 ) m 2 , 



while from ea. (|3.27|) . it can be inferred that 



= 7a(9 2 )A 



i f P(g 2 ) 

2g 2 



+ ia{ 9 2 ) A. 



(4.18) 



(4.19) 



(4.20) 



We notice the remarkable fact that the anomalous dimensions of the Gribov parameter and of 
the operator A 2 ^ are proportional to each other, to all orders of perturbation theory. 



It can now be checked that V is renormalization group invariant, namely 

ap 



(4.21) 



Finally, taking the derivatives of the action given in eq. (|4.14|) gives rise to 



1 dT 

A3 9A 



2 (iV 2 - 1) 3 (iV 2 - 1) 



g 2 N 



256vr 2 



m 



2 + Vm 4 - A 4 



Vm 4 - A 4 



In 



2 + Vm 4 - A 4 



'772 



2p 



+ 4 



m 



sf 



rrr 



A 4 



mr 



A 4 



m 2 - Vm 4 - A 4 8 
111 2^ + 3 



(4.22) 



'Using ln(z) = In \z\ + i arg(z) with — 7r < arg(z) < ty. 
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and 

dr 

dm 2 

Com 2 ( x _ Ci g2 \ + 3 (iV 2 - 1) 



9 2 \ Co / 256^2 



m 2 + Vm 4 - A 4 1 + , In 
/ V Vm 4 -A 4 7 



27i 2 



/ 9 / — : — \ / m 2 \ , m 2 — yJm A — A 4 o , 

+ 2 f m 2 - Vm 4 - A 4 J 11 , ] In — -m 2 



s/m*-\y 2Jl 2 3 



(4.23) 



4.2 Solving the gap equations. 



We have now all the ingredients at hand to search for estimates of the mass parameter m 2 and 
Gribov parameter A as solutions of the gap equations (|4.22|) and (|4.23j) . To avoid misinterpre- 
tations due to the suggestive use of the notation m 2 , we remark that, due to the presence of A, 
the mass parameter does not even appear as a pole in the tree level gluon propagator, see eq. (|6.2|) . 



Let us first consider the pure Gribov- Zwanziger case, i.e. we put m 2 = in the expression 
(|4.14|) . The relevant gap equation (horizon condition) reads 

OT 3 f 2(N 2 -1) 3 (N 2 - l) A 4 5\ 3(iV 2 -l) 



We remind here that the solution A = must be rejected. The natural choice for the renormal- 

T 

9 2 N 



ization scale is to set /i 2 = 4- to kill the logarithms, and we find 



16vr 2 



2_ A 2 
2 



4 . (4.25) 



In principle, from 



<?V) = ^s-, wi * h /3o = ^t^2, (4.26) 

MS 

ea. (|4.25[) could be used to determine an estimate for the Gribov parameter, however it might be 
clear that this is meaningless since the corresponding expansion parameter (|4.25f) is far too big. 

It is interesting to notice that, in a general massless renormalization scheme, the one-loop action 
with m 2 = would read 

{ N 2 - 1) 3A 4 {N 2 - 1) / _A_ 4 
2g 2 N 264vr 2 V n 47I 4 



r = - \. 9 ,/ A 4 ' (ln-4+o) , (4.27) 



with a an arbitrary constant. The corresponding gap equation equals 

9Y J 2{N 2 -l) 3 (N 2 - 1) / A 4 \ 3 (iV 2 - l) \ _ 

d\~ X { JN 64^ V ln i^ + °J " 64vr 2 J' ' (428) 

Denoting by A* a solution of eq. (|4.28|) , for the vacuum energy corresponding to (|4.27|) one finds 

F ro. ) 3(iy2 " 1} A * urn 

E ™ ~ T ^*> ~ 64vr 2 -J ■ ( 4 - 29 ) 
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This expression is valid for all /i and for all a. The vacuum energy is thus always nonnegative 
at one-loop order in the original Gribov-Zwanziger model. 

The gap equation (|4,23|) with A = obviously has the solution already obtained in [3J HHj 
where the restriction to the Gribov region £1 was not taken into account. We recall the values 

g 2 N 36 

- 0.193, (4.30) 



16tt 2 187 

9 II 

= ei2 

_3_ 17 

16^ 

tt2 _ ™2 



m 2 = eBA^«(2.031A w )^, (4.31) 
£vac = -^ef A^«-0.323A|g , (4.32) 



which were obtained upon setting /i = m to kill the logarithms. 

We shall now show that, in the MS scheme, the gap equations (|4.22|) - (|4.23|l have no solution 

with m 2 > when the restriction to the horizon is implemented (i.e. when A / 0). To this 
purpose, we introduce for m 2 > the variable 

t = 4 ■ (4-33) 

Evidently, we should only consider t > 0. 

Dividing the gap equations l|4.22|) - (|4.23|) by m 2 , they can be rewritten as 5 

16tt 2 3 ( „, m 2 2 1 , t , \ „„, 

21n— - + -+ In -hit , (4.34) 



5 2 iV 8 ^ 27i 2 3 ( X + 



24 /16vr 2 \ 322 m 2 4 2 - 1 , t , , 

+ — = 41n^--- In +21nt. (4.35) 



and 



13 \g 2 Nj 39 27i 2 3 (l + vT^T) 

where use has been made of the explicit values of Co an d Ci > which can be found in [51 1121 12*5] 

9 iV 2 - 1 161 iV 2 - 1 

The eans. ()4.35|) - (|4.36|) can be combined to eliminate ln|=7y, yielding the following condition 

68 /16vr 2 \ 161 t , i 

+ — = ^^ln 2=F(t). (4.37) 



39Vs 2 iV 39 (l + VT^T)' 

It can be checked that F(t) is real- valued and negative for t > 0, thus the r.h.s. of ea. (|4.37|) is 
always negative. Since the l.h.s. of eq. ([4.37|) is necessarily positive for a meaningful result (i.e. 
g 2 > 0), there is no solution with m 2 > 0. As already mentioned, there are a priori also possible 
solutions with m 2 < 0. 

To investigate the existence of a solution with m 2 < 0, it might be instructive to look again 
at the gap equations (|4.22j) and (|4.23|) from another perspective. We recall that, if the horizon 
is not implemented, i.e. A = 0, the gap equation (|4.23|) has two solutions, a perturbative one 

5 We have already factored out m 2 or A 3 since these are non-zero in the present case. 
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corresponding to m 2 = (no condensation) and a non-perturbative one corresponding to the 
m 2 given in ea. Q4.31[) . 



If we momentarily consider A as a free, adjustable parameter of the theory, ea, ()4.23|) dictates 
how m 2 becomes a function of the parameter A. From the result at A = 0, we could expect that 
two branches of solutions would evolve, one starting from the perturbative and one from the 
non-perturbative value of m 2 at A = 0. When A = 0, the choice for the scale ~p is quite obvious 

2 

from the requirement that all the logarithms In rqj- are vanishing. However, when A / 0, we 

notice that there are two kinds of logarithms present, being In : 
We opt to set 



+y5f^ and ln m2 -vS^ 



277 



M 2 



+ V' 



A 4 



(4.38) 



This reduces to /j, 2 = m 2 if A = 0, while it allows for 6 m 2 < 0. This is possible if m 4 < A 4 , as 
it was mentioned below ea. (|4.8|) . In this case, the size of both logarithms, In ™ 2 +^™?-^ an d 



In— — xi , is determined by their arguments, which are complex conjugate. 



2JZ A 



Let us specify to the case N = 3. In Figure 1, we have plotted the behaviour of m 2 {\ ). 
We see that next to the "non-perturbative" branch of solutions, starting from m 2 ^ 0, also a 
"perturbative" branch of solutions with m 2 < is emerging from m 2 = 0, in correspondence 
with our expectation. 



m A 2 



5 



-5 



•10 



50 100 150 200 250 



lambda A 4 



Figure 1: m 2 as a function of A 4 , in units 



However, A 4 is not a free parameter of the theory. We should require that A 4 is such that 
the doublet (A 4 ,m 2 (A 4 )) is a solution of the gap equation (|4.22|) . i.e. the horizon condition. In 
Figure 2, we have plotted the value of the horizon condition equation, as a function of A 4 . It 
is clear that no solution with m 2 > exists as the horizon condition never becomes zero. Of 
course, this is in correspondence with the foregoing general proof that there is never such a 
solution, independently of the choice of ~fi. However, we see that there is a single solution with 
m 2 < 0. 



'Evidently, /i should be real and positive, hence the modulus in eq.(|4 



16 



horizon condition 




Figure 2: The horizon condition (|4.22j) as a function of A 4 , in units A^g = 1. The top curve 
corresponds to the solutions of (|4.23|) with m 2 < and the lower curve to the solutions with 
m 2 > 0. 



The corresponding values for the expansion parameter, for the Gribov and mass parameter, 
as well as for the vacuum energy are found to be 

9 2 N 



16vr 2 



1.18 , (4.39) 



A 4 6.35lAlg, (4.40) 
m 2 -0.950A|3, (4.41) 
£ vac « 0.043Alg, (4.42) 



4.3 Intermediate comments. 



Although the MS expansion parameter (|4.39|) is too large to speak about reliable results, we 
nevertheless would like to raise some questions. Apparently, the solution of the coupled gap 
equations is laying on the "perturbative" branch, being the one with m 2 < 0. This gives rise 
to a positive value for the mass dimension two gluon condensate {Afy. When the restriction 
on the domain of integration in the path integral is not implemented, as in the previous papers 
[31 El EH] , {AV) was necessarily negative, the reason being that the action should be real- valued, 
as it was explained below ea. l|4.8jl . As already explained in the Introduction, a finding a bit 
unfortunate is that the vacuum energy is positive, ea. (14.4211 . which leads to a negative estimate 
for the gluon condensate <^„) via the trace anoma^CS- Essentially, we are thus left 
with the following questions: 

(i.) What is the sign and value of m 2 and thus of {A 2 ^? 

(ii.) What is the sign and value of E VSLC and the corresponding value for (^£^iF 2 J^l 

(hi.) Are these values better or not when the operator A 2 ^ is added to the original Gribov- 
Zwanziger model? 
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5 Changing and reducing the dependence on the renormaliza- 
tion scheme. 



We have already shown that the vacuum energy obtained in a one-loop approximation is always 
positive when the condensation of the operator A 2 ^ is left out of the discussion, using whatever 
renormalization scheme. 

To answer the foregoing questions (i.)-(iii.), one could investigate what happens at two- loop 
order. However, due to the already quite complicated structure of the one-loop effective action 
and to the fact that the calculations at higher loop order will not get any easier, this task is 
beyond the scope of the present article. Here, we shall mainly focus on the effects of a change of 
the renormalization scheme at the one-loop order. It could happen that, in a scheme different 
from the MS one, the vacuum energy is negative and/or that the coupling constant is small 
enough to speak about trustworthy results, at least qualitatively. 

Since to obtain an optimization of the renormalization scheme and scale dependence is a rather 
lengthy task, we shall not dwell upon technicalities in this section. The interested reader can find 
all details in Appendix B. We shall thus focus on the main results obtained after the optimization. 

Essentially, what we have done is replacing in the effective action (|4.14|) the quantities m 2 
and A 4 by their order by order renormalization scale and scheme invariant counterparts fh 2 and 
A 4 . The residual freedom in the choice of renormalization scheme can then be reduced to a single 
parameter 60 > related to coupling constant renormalization. As the vacuum energy is a physical 
quantity, it should in principle not depend on 60 . At the same time, the quantities fh 2 and A 4 
should be 60 independent by construction. This provides one with the interesting opportunity 
to fix the redundant parameter 6q by demanding a minimal dependence on it. 



The final one-loop action turns out to be given by 
r = 



Ar2_l) x- 26 A 4 ( X + B + (l-2b) f^ln^-60 



2N 



A A) 



+ «to V* ( X + A - £ + (1 - 2a) In * - bo ) ) + 3 ^ x 
2 V Co V ^ V/3o A) J J 256vr2 



fh 2 x~ a + V m 4 x" 2a - \ A x~ 2b 



^2 -a 

m x 



m A x 2a — X A x 26 




rh 2 x a + V fh A x~ 2a — \ A x~ 2b 5 



'Tr2 



277 



6 



m 2 x a — \f' 



m i x -2a _ \4 x -2b 5 

6 



'Tr2 



2JI 



(5.1) 



while the corresponding gap equations read 

1 dT 

A3 9A 

1 sr 



m 



2 Qfffi 



(5.2) 
(5.3) 



The definitions of all quantities appearing in the above expressions can be found in the Appendix 
B. 



In Figure 3, we collected the solutions of the scale invariant quantities fh 2 and A 4 as a function 
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of bo, while Figure 4 displays the vacuum energy E vac and the relevant expansion parameter, 
given by y = ■ For completeness, we have also shown the solutions which correspond to 
higher values of T and are as such describing unstable solutions. These are indicated with the 
thinner lines. 



I i 

m A 2 1*4 




Figure 3: The quantities m 2 and A 4 as a function of bo, in units = 1. 




Figure 4: The vacuum energy E vac and the expansion parameter y as a function of bo, in units 
A-ms = 1- 

5.1 Interpretation of the results. 

Let us first have a look at the plot of vacuum energy, on the l.h.s. of Figure 4. We notice 
that for bo < —0.33564...., the vacuum energy becomes negative. However, we cannot attach 
any definitive meaning to this result. In fact, as it can be seen from the Figures 3 and 4, the 
values of the vacuum energy and the supposedly minimally &o-dependent quantities in 2 and A 4 
are extremely &o-dependent. Very small variations in &o induce large fluctuations on e.g. the 
energy. This is indicative of the fact that the equations we have solved are not yet stable against 
bo- variations in the range of the values obtained for bo- The behaviour is better for, let us say 
bo > —0.2. However, in this case, we find again that the vacuum energy is positive. The vacuum 
energy E vac , as well as fh 2 and A 4 fall of to zero for growing bo- 

As an example, we set bo = 0, which corresponds to use the MS coupling constant. Then 
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we find, with the optimized expansion, 



y 


N 

= — — « 0.796, 
lbn z x 


(5.4) 


X 4 x~ 2b 


~ 7.939Aj^jg , 


(5.5) 


m 2 x~ a 


~ -O.8I4AI3, 


(5.6) 




« O.O63AI3, 


(5.7) 



results which are in fair agreement with the naive MS results (|4.39j l- (|4.42jl . We notice that the 
expansion parameter y is already smaller than 1, but still relatively large, while the vacuum 
energy is indeed positive. 

The conclusion than can be drawn from this section is that we cannot find a reliable result 
with negative vacuum energy and hence positive gluon condensate (F 2 U ) using a one-loop ap- 
proximation. We see therefore that, in order to be able to give a reasonable answer to the 
questions concerning the sign of m 2 and E vac and to get more trustworthy numerical values, the 
two-loop evaluation of the effective action T, at least in the MS scheme, would be very useful. 

6 Consequences of a non-vanishing Gribov parameter. 

Before turning to the final conclusions, we shall give in this section a brief account of some well 
known consequences stemming from the presence of the Gribov parameter, to emphasize the 
important role of this parameter. 

6.1 The gluon propagator. 

If there is no generation of a mass parameter due to \Afy, we can consider just the action (|2.2[) . 
Then the tree level gluon propagator turns out to be 

M,-**!^ ■ ((U) 

This result, first pointed out in |2B], was obtained by retaining only the first term of the nonlocal 
horizon function (|1.4|) . corresponding to the approximation — dD ~ — d 2 . The gluon propagator, 
eq. ([6.1|) . is suppressed in the infrared region due to the presence of the Gribov parameter A. In 
particular, the presence of this parameters implies that (A^Afy vanishes at zero momentum, 
p = 0. When the possibility of the existence of a dynamical mass parameter in the gluon 
propagator is included, by investigating the condensation of A 2 ^ the tree level gluon propagator 
reads 

/j.jt) s ( ,„ _ a*) = ^ f >4 ( J( ,„ _ a*) . (6 . 2 ) 

\ M / p p 2 V P J p 4 + m 2 p 2 + V P ) 

This type of propagator is sometimes called the Stingl propagator, from the author who used it 
as an anzatz for solving the Schwinger-Dyson equations, see |66| for more details . 

However, it should be realized that ea. (|6.2|) describes only the tree level gluon propagator. 
In particular, to produce a plot of the form factor D(p 2 ) as a function of the momentum p, 
which would allow to make a comparison with the results obtained in lattice simulations, see 
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e.g. 31 for N = 3 and [321 OH for N = 2, one should go beyond the zeroth order approxi- 
mation, for example by including higher order polarization effects and/or trying to perform a 
renormalization group improvement. In general, these corrections will also be dependent on the 
external momentum p. 

6.2 The ghost propagator. 

Even more prominent is the influence of the Gribov parameter on the infrared behaviour of the 
ghost propagator, which can be calculated at one-loop order using the modified gluon propagator 
(|6.1|) or 1)6.2(1 with their respective gap equations (|1.6|) and (|4.3jl . In both cases, the infrared 
behaviour of the ghost propagator [JBl E3 E3 03 E2j is shown to be 



The original Gribov- Zwanziger model corresponds to m 2 = 0. Thus, the ghost propagator is 
strongly enhanced in the infrared region compared to the perturbative behaviour, if the re- 
striction to the first Gribov region is taken into account. It is important to notice that this 
behaviour of the ghost propagator is preserved in the present treatment, due to the peculiar 
form of the gap equation (|4.3j) implementing the horizon condition. In particular, from the 
expression for the effective action in ea. (|4.2|) . one sees that, while the term quadratic in the 

2 

field a, i.e. 2g^' contains the LCO parameter the first term which depends on the Gribov 
parameter, i.e. — 4(N 2 — 1)7 4 , does not contain any such new LCO parameter. This important 
feature follows from the fact that no new parameter has to be introduced in order to renormalize 
the term (M^Vf^ 1 — U® l N^, as remarked in ea. (|A.35|) . While the parameter £ is required to 
take into account the ultraviolet divergences of the vacuum correlator (A^ 1 (x)A 2 (y)) , which are 
proportional to r 2 , no such a parameter is needed for [M^V™ — UjfN™) which, upon setting 
the external sources to their physical values, gives rise to term — 4(iV 2 — 1)7 4 in the expression 
(|4.2|) . Said otherwise, this term is not affected by the presence of a new parameter which would 
be required if ea. (|A.35|) would not hold. As a consequence, the factor "1" appearing in the 
left hand side of the gap equation (|4.3|) is, so to speak, left unchanged by the quantum correc- 
tions. It is precisely that property which ensures, through a delicate cancelation mechanism, see 
[13 EH 03 EH > the infrared enhancement of the ghost propagator. 

Analogously to the case of the gluon propagator, a more detailed study of higher order cor- 
rections would be needed in order to obtain a plot of the ghost form factor G{p 2 )- 

6.3 The strong coupling constant. 

Usually, a nonperturbative definition of the renormalized strong coupling constant an can be 
written down from the knowledge of the gluon and ghost propagators as, see e.g. |39l I34| 



where T> and Q stand for the gluon and ghost form factors as defined before. This definition 
represents a kind of nonperturbative extension of the perturbative results ()A.33|) . According to 




(6.3) 



where J stands for the real, finite integral given by 




(6.4) 



oir{p 2 ) = a R (v)V(p 2 ,LL)g 2 (p 2 ,{i) 



(6.5) 



21 



Schwinger-Dyson studies E21 ESI EJ EI3 EE] > those form factors satisfy a power law behaviour 
in the infrared 

limP(p 2 ) oc (p 2 ) 9 , 

limg(p 2 ) oc (p 2 Y , (6.6) 

where the infrared exponents 6 and u obey the sum rule 

9 + 2ui = . (6.7) 

Such a sum rule suggests the development of an infrared fixed point for the renormalized coupling 
constant, ()6.5|) . as also pointed out by lattice simulations for the SU{2) as well as for the SU(3) 

case jMHSHEEl, 

lim a(p 2 ) = a c . (6-8) 

The existence of a fixed point in this reasoning is dependent on the sum rule rather than on 
the precise value of the exponents. We refer to the already quoted literature for more details on 
the value of these exponents. We end by noticing that the form factors of the gluon and ghost 
propagator in our zeroth order approximation give rise to the sum rule (|6.7|) , since we have 9 = 2 
and oj = —1. Moreover, without Gribov parameter, the sum rule (|6.7|) is lost, and thus there is 
no indication for an infrared fixed point. 



6.4 Positivity violation. 

The behaviour of the gluon propagator is sometimes used as an indication of confinement of 
gluons by means of the so called positivity violation, see e.g. |671 168j and references therein. 

Briefly, when the Euclidean gluon propagator D(p) = V ^ - is written through a spectral repre- 
sentation as 

the spectral density p(M 2 ) should be positive in order to have a Kallen-Lehmann representation, 
making possible the interpretation of the fields in term of stable particles. We refer to |67| I68j 
for more details. One can define the temporal correlator 



r+oo 

C(t) = / dMp(M 2 )e- Mt , (6.10) 
J o 

which is certainly positive for positive p(M 2 ). The inverse is not necessarily true. C(t) can be 
also positive for a p(M 2 ) attaining negative values. However, if C(t) becomes negative for certain 
t, then a fortiori p(M 2 ) cannot be always positive. Using a contour integration argument, it is 
not difficult to show that C(t) can be rewritten as 

i r+oo 

C{t) = — / e-^D(p)dp . (6.11) 
Let us consider the function C(t) using the tree level propagator l|6.2j) . thus using 

2 

D(p) = 1 A4 • (6.12) 

+ p L m L + — 

We can consider several cases 7 : 

7 Each of the following expressions for C(t) is obtainable via contour integration. 
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if A = (thus m > 0), one shall find that 



-mt 

C(t) = — . (6.13) 
2m 



This function is always positive, 
if m 2 = 0, 



Lt 



e 2 / Lt Lt\ , 
C(t) = — - cos sin— , (6.14) 



2L \ 2 2 J 

and clearly, this function will attain negative values for certain t 

in any other case, the correlator C(t) is found to be 

C (t) = I [ ^ e -yfint + e -V^ 
2 \_U)\ — UJ2 UJ2~ 0J\ 

where the decomposition 

p 2 L0\ 1 U>2 1 



(6.15) 



(6.16) 



P 4 + p 2 m 2 + ^ - W2 P 2 + 0J\ Ul — L02 P 2 + L02 

has been employed. It is understood that y/uJ\ (y^i) is the root having a positive real part. 



In 

'2_ 



If we assume that m 4 > A 4 , then u)\ > 0J2 and C(t) becomes negative for t > ,, 
In the case that m 4 = A 4 , or u;i = 0^2, one finds that C(t) = e 4 . ^ (1 — y / tiJit), which 



can also become negative. If m 4 < A 4 , we can reintroduce the complex polar coordinates 
R and for the complex conjugate quantities uj\ and 002- If cos f > 0, ea. (|6.15|) can be 
rewritten as 

By choosing an appropriate value of t > 0, also this expression can be made negative. An 
analogous expression and conclusion can be derived in case that cos | < 

We conclude that, when the restriction to the Gribov region 0, is implemented, the function C(t) 
exhibits a violation of positivity when the tree level propagator is used, with our without the 
inclusion of (^4 2 t )- 

The goal of this section was merely to provide some interesting consequences when the restric- 
tion to the first Gribov region Q is implemented. Higher loop effects, which shall be momentum 
dependent, would also influence the behaviour of the gluon and ghost propagator. Hence, to 
give a sensible interpretation of the behaviour of e.g. the form factors and of the strong coupling 
constant an, a more detailed analysis than a tree level one is necessary. This is however far 
beyond the aim of this work. 



7 Conclusion. 

In this work we have considered SU(N) Euclidean Yang-Mills theories in the Landau gauge, 
d^A^ = 0. We have studied the condensation of the dimension two composite operator A 2 when 
the restriction to the Gribov region f2 is taken into account. Such a restriction is needed due to 
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the presence of the Gribov copies [IE], which are known to affect the Landau gauge. 

In a previous work [inj> the consequences of the restriction to the region f2 in the presence 
of a dynamical mass parameter due to the gluon condensate (A^Afy were studied by following 
Gribov's seminal work |3H1- Here, we have relied on Zwanziger's action |51[I52|. which allows to 
implement the restriction to the Gribov region Vl within a local and renormalizable framework. 
We have been able to show that Zwanziger's action remains renormalizable to all orders of per- 
turbation theory in the presence of the operator A^, introduced through the local composite 
operator technique [HI El 1201 ESI • The effective action, constructed via the local composite 
operator formalism [H] obeys a homogeneous renormalization group. The explicit form of the 
one-loop effective action has been worked out. We have seen that, considering the original 
Gribov-Zwanziger model, i.e. without including the operator A^, the vacuum energy is always 
positive at one-loop order, independently from the choice of the renormalization scheme. A 
positive vacuum energy would give rise to a negative value for the gauge invariant gluon con- 
densate (i* 1 /^,), through the trace anomaly. Furthermore, by adding the operator A^, we have 
proven that there is no solution of the two coupled gap equations at the one-loop order in the 
MS scheme with (Afy < 0. Nevertheless, when {A 2 ^ > 0, a solution of the gap equations was 
found, although the corresponding expansion parameter was too large and the vacuum energy 
still positive. 

In order to find out what happens in other schemes, we performed a detailed study, at low- 
est order, of the influence of the renormalization scheme. We have been able to reduce the 
freedom of the choice of the renormalization scheme to two parameters, namely the renormal- 
ization scale ~p, and a parameter 6o, associated to the coupling constant renormalization. We 
reexpressed the effective action in terms of the mass parameter fh and Gribov parameter A, 
which are renormalization scheme and scale independent order by order. The resulting gap 
equations for these parameters have been solved numerically. Although a solution with negative 
vacuum energy was found, we have been unable to attach any definitive meaning to it. This is 
due to the fact that the results obtained turn out to be strongly dependent from the parameter 
&o- This brought us to the conclusion that we should extend our calculations to a higher order 
to obtain more sensible numerical estimates. 

The mass parameters fh and A are of a nonperturbative nature and appear in the gluon and 
ghost propagator. Even if we lack reliable estimates for these parameters, some already known 
interesting features can be recovered. For a nonzero mass and Gribov parameter, there is a qual- 
itative agreement with the behaviour found in lattice simulations and Schwinger-Dyson studies: 
a suppressed gluon and enhanced ghost propagator in the infrared, while further consequences 
of the Gribov parameter are e.g. the possible existence of an infrared fixed point for the strong 
coupling constant and the violation of positivity related to the gluon propagator. 
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A Appendix A. 



In this Appendix, we have collected all details of the multiplicative renormalization of the 
Zwanziger action in the presence of the operator A 2 . 



A.l Ward identities. 



In order to begin with the algebraic characterization of the most general counterterm needed 

for the renormalizability of the complete action E of ea. ()3.1[) . let us first give the set of Ward 
identities which are fulfilled by E. These are 

• the Slavnov- Taylor identity 

5(E) = 0, (A.l) 

with 

, . f l4 / <5E SE <5E <5E , 5E „ <5E „ 5E 

• the Landau gauge condition and the antighost equation 

H = MS. ( A - 3 ) 



5E 5E 



+ 9,j^ = 0, (A.4) 



• the ghost Ward identity 



S a E = AS , (A.5) 



with 



(A.6) 



and 



Kl = 9 J d A xf abc - L b c c ) . (A.7) 

Notice that the term A^, being linear in the quantum fields A®, c a , is a classical breaking, 
the linearly broken local constraints 

ft 

fj, fl 



25 



the integrated Ward identity 

faff^+^te+U*™) (A. 12) 

V Suj m 5c a 5Kf, ' 

the exact IZij symmetry 

KijE = , (A. 13) 

with 

K « - 1 d ' x - + v "sk - "'"^ ■ (A - 14) 



A. 2 Algebraic characterization of the counterterm. 

Having established all the Ward identities fulfilled by the complete action E, we can now turn 
to the characterization of the most general allowed counterterm S c . Following the algebraic 
renormalization procedure [2^, E c is an integrated local polynomial in the fields and sources 
with dimension bounded by four, with vanishing ghost number and Qj-charge, obeying the 
following constraints 



8X C £S C 

+ ^Jn™ ~ ' 

8E „ J£ 

+ d u — — = 0, 



8Tp ai ^8M a 



<5£ c 



and 



(A.15) 



g a S c = 0, (A.16) 

/ ^ / c a^_ +lJ ai^ U?^-) = , (A.17) 

TZijT, = , (A. 18) 



23 s £ c = , (A.19) 
where By, is the nilpotent linearized Slavnov- Taylor operator 

+ *sr + + M;i w + + T i) • (A - 20) 
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£s£ s = . (A.21) 

As it was shown in |511 1521 153] . the constraints (|A.15|) imply that S c does not depend on the 
Lagrange multiplier b a , and that the antighost c a and the i-valued fields ipf, ujf, Tp^, uJ" can enter 
only through the combinations 

K = K- + d^- g r bc u b ;^-gf ahc v b ^\ 
v? = C + V. 

M« = V? + . (A.22) 

Therefore, S c can be parametrized as follows 

S c = S C (A) + J (fx (a ig f abc L a c b c c + a 2 K^c a + a 3 g f abc K^c + a 4 f abc V™U^c c 

+ a 5 V™M™ + a 6 U™N™ + ^TAZA« + ^t 2 + a^A^c* + a m c a dA a ) , (A.23) 

where S C (A) depends only on the gauge field A a ^, and with a\, a±o arbitrary parameters. 

Notice, however, that there is no mixing in expression (JA.23JI between MT, N™, V™, and 
the sources r, r\. This is due to the dimensionality and to the Qy-charge. It is precisely the 
absence of this mixing that will ensure the renormalizability of the Zwanziger action in the 
presence of the composite operator A®A^. From the ghost equation (|A.16|) it follows 

ai = a 3 = oio = , 

04 = -g{a 6 + a 5 ). (A.24) 
^From the equations (|A.17|) and (|A. 18|) we obtain 

a 6 = -a 2 . (A.25) 

Finally, from ea. (|A.19|) it turns out that 

a 5 = a 2 , 

a g = a 7 - a 2 , (A. 26) 

and 

S C (A) = a S YM +a 2 J d^ 6 -^ . ( A .27) 

In summary, the most general local invariant counterterm compatible with all Ward identities 
contains four arbitrary parameters, do, a 2 , 07, as, and reads 

S c = a S YM + a 2 f d*x + K^c* + V«M« - U™N™ 

+ / d 4 x r^A« + ^Cr 2 + (or - a 2 ) n A%d^) . (A.28) 



27 



A. 3 Stability and renormalization constants. 

Having determined the most general local invariant counterterm S c compatible with all Ward 
identities, it remains to check that the starting action S is stable, i.e. that S c can be reabsorbed 
through the renormalization of the parameters, fields and sources of S. According to expression 
(|A.28|) . S c contains four arbitrary parameters ao, a 2 , 07, a§, which correspond in fact to a 
multiplicative renormalization of the gauge coupling constant g, the parameters £, and of the 
fields (f> = (A*, c a , c a , b a , tpf, ojf, Tpf, uf) and sources $ = {K a ^, L a , M™, Nf, Vj*, Uj*,r, rj), 
according to 

S( ff ,C,0,f) + ??S c = S( 5o ,Co,0o,$o) + O(r/ 2 ), (A.29) 

with 

9o = Z g g , Co = Z ( C , (A.30) 

and 

3> = . (A.31) 

The coefficients ao , a 2 are easily seen to be related to the renormalization of the gauge coupling 
constant g and of the gauge field A a , 

z, = (1 + 1?) ■ 

Z'f = (l + „(a 2 -f)) . (A.32) 

^From expression (|A.28|) it follows that the Faddeev- Popov ghosts (c a ,c a ) and the i- valued fields 
((fi,Luf,Tpf ,ZJf) have a common renormalization constant, given by 

Z c = Zc = Z v = Z ¥ = Z u = = (1 - 7]a 2 ) = Z g 1 Z A 1/2 . (A.33) 

Ea. ()A.33|) expresses a well-known renormalization property of the Faddeev-Popov ghosts (c°,c a ) 

in the Landau gauge, stemming from the transversality of the gauge propagator and from the 

factorization of the ghost momentum in the ghost-antighost-gluon vertex. We see therefore that, 

in the present case, this property holds for the i-valued fields ((pf,uf,lp^,LJ^) as well. Similarly 

to the ghost and the i-valued fields, the renormalization of the sources (M® 1 , N™, V^*, Ujf) is 

1/2 

also determined by the renormalization constants Z g and Z A , being given by 

Z M = Z N = Z v = Z v = Z~y 2 Z- 1/A . (A.34) 

It is worth noticing here that equation (|A.34|) ensures that the counterterm a 2 (V^ l Mff — U™N™~) 
can be automatically reabsorbed by the term (— M^ l V™ + U™N™^ in the expression H3.5|) with- 
out the need of introducing new free parameters. Indeed, 

-M V = -MVZ 2 M = -MVZ g l Z A 1 ' 2 = -MV + ea 2 MV . (A.35) 

Concerning now the parameters aj, as, they are easily seen to correspond to a multiplicative 
renormalization of the local source r and of the parameter £, according to 

t = Z t t , Z T = 1 + 77(07 - 2a 2 + ao) , 

Co = Z c ( , Z c = l + 7](-a 8 - 2a 7 + 4a 2 - 2a ) . (A.36) 
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Moreover, we would like to underline that there exists even an extra relation, namely 

Z T = Z g Z A 1/2 . (A.37) 

It can be proven by introducing the operator A 2 ^ through a more sophisticated set of local 
sources, like it was done in We will not repeat that analysis here, we only mention that a 
key ingredient in the proof of relation ( IA.37|) was the presence of the ghost Ward identity, and 
since the Zwanziger action possesses that identity, eq. (|A.5|) . one can proceed along the lines of 
jllj . Thus, there are in fact only three independent renormalization factors present. 

B Appendix B. 

In this Appendix, we give the detailed analysis of the procedure used to optimize the renormal- 
ization scheme and scale dependence, which was summarized in section 5. 

B.l Preliminaries. 

Before coming to the actual computations, let us first discuss some results which will turn out 
to be useful. 

Consider again the action S of eq. (|3.8|) . Due to the rich symmetry structure of the model, 
encoded in the Ward identities (|A. 1 j) - (|A. 141) . and due to the extra relation (|A.37|) . only three 
renormalization factors remain to be fixed, namely Z g , Za and Z^. Apparently, this means 
that we would need three renormalization conditions in order to fix a particular renormalization 
scheme. However, taking a look at the bare action associated with expression eq. (|3.8j) . we would 
find the following relations 

Co = z c (, 

CoT D 2 = W'Z^T 2 , 

t = Z t t , (B.38) 

from which it follows that 

Z c ( = PCX . (B.39) 

Since the bare quantity Q is renormalization scheme and scale independent and since £ always 
appears in the combination Z^ in the action, it follows that only Z g and Za are relevant for 
the effective action, because Z T can be expressed in terms of these two factors. Consequently, 
we would only need two renormalization conditions to fix the scheme. Obviously, we can equally 
well choose to make use of, for example, Z g and Z T as the two independent renormalization 
factors, corresponding to coupling constant and mass renormalization. 

We will change from the MS to another massless renormalization scheme by means of the fol- 
lowing transformations 8 

g 2 = g 2 (i + b g 2 + b ig 4 + ---) , 
A = A(l + c 5 2 + ci5 4 + ---) , 
m 2 = m 2 (1 + d g 2 + di^ 4 + • • • ) , (B.40) 
8 Barred quantities refer to the MS scheme. 
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where the parameters hi, q and di label the new scheme. However, we should keep in mind 
that the renormalization of the Gribov parameter A is not independent of that of g 2 and m 2 . 
Eliminating 7a (ff 2 ) between eqns. (|4.19[) and (|4.2()j) . yields 



i\(g 2 



(B.41) 



This relation, valid to all orders of perturbation theory, implies the existence of relationships 
between the coefficients b{, Ci and di. For further use, we shall explicitly construct the relation 
between bo , cq and do. Let us adopt as parametrization of f3(g 2 ), 7 m 2 (g 2 ) and 7a (<? 2 ) 

(3(g 2 ) = -2(/3 9 4 + /?i9 6 + ---) , 



i\(g 2 



log 2 + 7i5 4 + • 
AoS 2 + Ai/ + 



(B.42) 



and an analogous one in the case of the MS scheme. Then, one computes 
*W = /^[A(lW + ..•)] 



= A (\ g 2 + (Ai + c A - 2/3 c ) g 4 + • • • ) , 
which can be expressed in terms of 7, and /3j by exploiting the relation (|B.41|I . We find 



dX 
cTp, 



-2/3 ~ 70 2 , 
3 5 + 



-2&-71 , -2A) - 70 » 4^ 
: hc 2/3 c 5 + 



(B.43) 



(B.44) 



We can also calculate /x^= by first exploiting the relation (|B,41|) . obtaining 



dA 
5/1 



[(-2/? - 7o)5 2 + (-2A - 7i)5 4 + • • • ] [A (l + c g 2 + •••)] 

[(-2/? - 70) <? 2 + (co(-2/3 - 70) - 2#l - 71 - 2/3 (-d + &o) )/ + ••• ] (B.45) 



In the previous expression, we had to express 7^ in terms of 71; a task accomplished by using 
the relation 

7i = 71 - 2(3 d - 7060 , (B.46) 

which can be obtained along the same lines of the previous calculations. It should also be noted 
that 70, Po and (3\ are renormalization scheme independent quantities. Thus, the identification 
of eans. ()B.44|) and ()B.45|) gives the desired relation, given by 

1 



co = t(&o - d ) ■ 



(B.47) 



We now perform the transformations (|B,4fl|) on the action (|4.14|) . which was calculated in the 
MS scheme, to obtain it in a general scheme. 



r 



+ 



+ 



(N 2 - 1) A 4 
2g 2 N 
3 (N 2 - 1) 
256vr 2 

(m 2 — \/ 



(1 + 4c o5 2 - b g 2 ) + 



Com 4 



1 



A 4 



hi 



2g 2 

m 2 + Vm 4 - A 4 



^g 2 + 2dog 2 - bog 2 
Co 



>tt2 



2/7 



5 



777 ' 



A 4 



In 



771 



2/r 



1 -A 4 _ 5 
6 



(B.48) 
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while the gap equations now read 



8T 



2(iV 2 -l) „ 3(iV 2 -l)A 3 

; -A 3 (l + 4c g 2 -b g 2 + 1 



9 2 N 

( 2 



256?r 2 



Vm 4 - A 4 



4 -A 4 | ,,.,2 
In 



2JL 2 



4 _ _\4 i /J' - 

+ 4- 



Vm 4 - A 4 



ire* — A 4 



ki- 



rn 



A 4 



! 7?2 



2/7 



or 

dm 2 



5 2 V Co 



+ 



3 iV 2 



1 



256vr 2 
2 fm 2 - 



^m 2 + x/m 4 



A 4 



m 



Vm 4 - A 4 



In- 



m 



2 + _ A 4 



27I 2 



A 4 



mr 



A 4 



In 



Vm 4 



A 4 8 



tt2 



2/7 



-m 
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We mention that, in the case in which m 2 > 0, similar algebraic manipulations as those leading 
to the condition (|4.37|) . give a more general equation 



68 /16vr 2 \ 161 16vr 2 /32 68 
39 \gW ) + ~W \Y C ° ~ 39 

or, using the relation (|B.47|) . 

68 /16vr 2 \ 161 16vr 2 /12 
39 J + lS9 + ~N~ V 13 ° 
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:hl 



(i + vT^I) 



2 ' 



(B.50) 



176 
~39~ ( 



In 



VT^t (i + t) 



2 ' 



(B.51) 



^From this expression, it is apparent that a sensible solution with m 2 > might exist, depending 
on the values of the renormalization parameters do (~ mass renormalization) and &o (~ coupling 
constant renormalization) . 

Frequently used are the so-called physical renormalization schemes whereby, loosely speaking, 
one demands that the quantum corrected quantities reduce to the tree level values at a certain 
scale 71. However, it turns out that such an approach is not particularly useful to implement 
in the current case due to the presence of the several scales. Therefore, the question arises 
how one can make a somewhat motivated choice for the arbitrary parameters, labeling a certain 
renormalization scheme. In the next subsection we shall discuss a way to reduce the freedom in 
the choice of the renormalization parameters. The method relies on the possibility of performing 
an optimization of the renormalization scheme dependence, as illustrated in |70[ I71j. 

B.2 Optimization of the renormalization scheme. 



Consider a quantity g that runs according to 

_dg 



where 



^=7,(5 )Q, 



(B.52) 
(B.53) 
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To g, we can associate a quantity g that does not depend on the choice of the renormalization 
scheme and which is scale independent. It is defined as 



(B.54) 



whereby 



(B.55) 



It is apparent that g will not depend on the scale 71. It can also be checked |7U1 171| that g is left 
unmodified by a change of the renormalization scheme, implemented through transformations 
analogous to those of eqns. (IB~40l) . The equation (|B~55l) can be solved in a series expansion in 
g 2 by noticing that 

- ' ' ' %F e (9 2 )- (B.56) 



Tlj=F g (g 2 ) = P(g 2 )— T e (g 2 ). 



dg 2 

Then, the above differential equation can be solved in a series expansion in g 2 , more precisely 
by 



T e ,0 



W) = (ff ) 2/30 1 + 



1 (lg,l PxIq, 



v+- 



(B.57) 



2 V Po pi 

Consider once more the MS action T given in eq. (|4.14|) . We shall now replace the MS variables 
m 2 and A by their renormalization scheme and scale independent counterparts m 2 and A, which 
are obtained as before. By inverting ea. (|B.57|) . one has 



m 



Moreover, introducing the notations 



A 



7o 

2/?o ' 
7i _ ffiTo 
Po Pi 



(sT^(i-~(^ 



2 v/3 /? 2 r 
-if) 9 + 



the one-loop action is rewritten as 



2 V A 



t A 

6 = _ /Sb' 



m 



A . 



5 



Po Pi 



(B.58) 
(B.59) 

(B.60) 
(B.61) 



2N yy ' \ g 2 J 2 



ff 2+A Co/ 2567T2 



In- 



m 2 (5 2 ) a + JmH~9 2 ) 2a - A 4 (g 2 ) 2b 



Tr2 



2/7 



+ ( m 2 (5 2 ) a - Jm*(g 2 ) 2a - \ A {g 2 ) 2b I I In ■ 



m 2 {g 2 ) a - \/m*(g 2 ) 2a - A 4 (s 2 ) 2b 

tt2 



277 



(B.62) 
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The action (|B,62|) is still written in terms of the MS coupling g 2 . Performing the first transfor- 
mation of (|B.40|) . T can be reexpressed as 

{N2 ~\ 2 ) 2b ^(^ + B-b + 2bb 



2N 



+ ^m\g 2 ) 2a [^ + A-b + 2ab - ^ ) + 



9 



Ci 
Co 



rn \g 2 ) a + dm\g 2 ) 2a - A 4 ( 5 2 ) 
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hi 



3 (N 2 - 1) 
256?r 2 

m 2 (g 2 ) a + \/m i (g 2 ) 2a - \ A {g 2 ) 2b 

~W 2 



+ 



fh 2 (g T - ^m\g 2 ) 2a - X 4 (g 2 ) 



2\2b 



111 



m 2 (g 2 ) a - ^m A {g 2 ) 2a - X 4 (g 2 ) 
2p 2 



2^26 



(B.63) 



So far, we have constructed an action which is written in terms of renormalization scale and 
scheme independent variables A and m 2 and the coupling constant g 2 (77). This is a certain 
improvement, since we are not faced anymore with a choice of the parameters di, related to the 
renormalization of the Gribov and mass parameter. The remaining freedom in the choice of the 
renormalization scheme resides in the coupling constant, labeled by the parameters bo,bi,..., 
and in the scale fl. Of course, the higher order coefficients bi, i = 1, . . . do not show up here, 
since we have restricted ourselves to the one-loop level. Nevertheless, we will perform one more 
step, since the dependence on the coupling constant renormalization can be reduced to solely 
bo, by expanding the perturbative series in inverse powers of 



x = Po ln J2 ' 



rather than in terms of g 2 . For another illustration of this, see e.g. |7D| I71j . 
constant g 2 can be replaced by x since g 2 is explicitly determined by 



(B.64) 



The coupling 



9 



0! ln & 




(B.65) 



In [22], the relation between the scale parameter A, corresponding to a certain coupling constant 
renormalization, and that of the MS scheme, A^g, was found to be 



b 

A = e 2 "oAt 



(B.66) 



One finally gets the expression (|5.1j) . We notice that this alternative expansion is correct up to 
order (^) . 

In principle, we can solve the two equations (|5.'2j) - (|5.3|) for the two quantities m* and A*, which 
will be functions of the two remaining parameters 71 and fro- However, by construction, we know 
that fh as well as A should be independent of the renormalization scale and scheme order by 
order. This gives us an interesting way to fix these parameters by demanding that the solutions 
m*(/Z, bo) and A* (77, &o) depend minimally on bo and 77. Since this would give a quite complicated 
set of equations to solve, we can make life somewhat easier by reasonably choosing the scale 9 77 



This can be motivated thanks to the scale independence of the ^-quantities. 
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in the gap equations (|5.2[) - l|5.3j) . In analogy to the choice for fi 2 done in the previous equation 
()4.38j) . we shall now set 



M 2 



m 2 x a + V m A x~ 2a — \ A x~ 2b 



(B.67) 



In order to proceed, we still have two quantities at our disposal to fix the remaining parameter 
&o- In fact, we can also take the vacuum energy E vac in consideration since, being a physical 
quantity, it should depend minimally on the renormalization scheme and scale. Therefore, we 
could determine the value for bo by demanding that 



T(6 



+ 



dfhi 



db 



+ 



dE v 



db n 



(B.68) 



is minimal w.r.t. the parameter bo. This seems to be a reasonable candidate. When its depen- 
dence on bo is small, then the dependence of m, A and £\,ac on bo is necessarily small too. The 
ideal situation would be that T is zero for a certain bo- If no such an ideal bo would exist, we 
weaken the condition by requiring that T is as small as possible. The condition ()B.68|) to fix 
bo can be considered as some kind of principle of minimal sensitivity a la Stevenson 72 . An 
alternative that is sometimes used is a fastest apparent convergence criterion, where it is de- 
manded that the quantum corrections are as small as possible compared to the tree level value. 
For example, if we denote by rM the action to order (^) and by rM to order (^) , we could 
demand that 



(B.69) 



is as small as possible when the parameters fulfill the gap equation describing the vacuum of 
the theory. 

Before continuing with explicit calculations, let us just remark here that the other logarithm, 



namely In • 



— , could become large for a small argument, thus when X A x 2b 
would be small compared to m A x~ 2a . However, it is harmless since it appears in the form of 



uhiu, while we know that ulnul 



0. 



B.3 Numerical results. 

Let us first give some numerical factors we need. From e.g. [HS], we infer that 

3 N 



01 



34 / N 
Y V 16vr 2 



7o 



2 16vr 2 ' 



71 



95 ( N V 
24 V 16^ 2 J 



and hence, from the relation (|B.41|) . 

35 N 

An 



A i 



449 / N V 
~96" V 16tt 2 J 



24 16tt 2 ' 

This means that, for any N, the quantities a and b in eq. ([B.60|) are found to be 
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44 



(,=1. 



(B.70) 



(B.71) 



(B.72) 
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It is instructive to consider once more the original Gribov-Zwanziger model by setting m = 
and by solving the gap equation (|5.2|) . If A* is a solution of this equation, then it is not difficult 
to show that the corresponding vacuum energy is given by 

3(iV 2 - 1) \j 

E ™ ~ 64^ T ' (B - 73) 

for any choice of ~p 2 . Thus, also with the improved perturbative expansion, the vacuum energy 
of the original Gribov-Zwanziger is always nonnegative at the lowest order. 

Let us return to the model we were investigating. We solved the gap equations stemming 
from (|5.2|) - (|5.3|) numerically. 

Let us first search for a possible solution of the gap equation in the region of space deter- 
mined by m 4 x _2a > \ A x~ 2b . Taking a look at the action ()5.1|) . it might be clear that the gap 
equations derived from it will be coupled and hence quite complicated to solve numerically. 
From the calculational point of view, it is useful to introduce new variables, defined by 



fn 2 x a + 



\/fh i x~ 2a - A 4 x" 26 



ux = , (B.74) 



m 2 x 



^2 = g ' (R75) 

with the inverse transformation 

m 2 x~ a = u)\ + UJ2 , 

A 4 x" 2fe = 4wiw 2 • (B.76) 

This defines a mapping from the space m 4 x _2a > A 4 x~ 2b > to uj\ > 0J2 > 0. One checks that 
the gap equations l|5.2|) - (|5.3jl are equivalent to 

T(ui,u)2) = 0, (B.77) 



Id Id 



U>\ — OJ2 du>i U>1 — U>2 dbJ2 



r(wi,wa) = 0. (B.78) 



We notice that the case in which uj\ and 002 would become equal, i.e. m A x~~ 2a = X 4 x~ 2b , should 
be treated with some extra care. Let us therefore first assume that u\ > u)%. Then the two 
equations ()B.77|) - I|B.78|) can be recombined to 

d 

T = 0, (B.79) 



d 



The action T(ti)\,u)2) is explicitly given by 



. (B.80) 



(N 2 -l) Co , 

-2 V - J U 1 u 1 uj2 + ^U 2 (wi + ^ 2 ) 2 



3 (N 2 - 1) 
64^ 



W~S) +UJ2 { W "6 



(B.81) 
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where 



V 1 = x + B + (l-2b)(^ln^--b ) , (B.82) 

VPi) Po / 

U2 = x + A-^ + (l-2a)(^\n^--b ) . (B.83) 
It is not difficult to work out the gap equations ()B.79|) - (|B.8fl|) . being given by 



iV ^ v ' 32vr 2 V 3 



-2^0^ + Co^ + o») + 3(iV 3 2 2 ; 2 lV2 (~~ + In g ) = 0. (B.85) 



/,From the explicit expression of the gap equations and of the action itself in terms of ui\ and 
LU2, the advantages of using these variables should be obvious, since we can decouple the two 
gap equations. Explicitly, since 71 2 =(*)%, one finds from ea. (|B.84|) . 



w 2 = fgh, oJi , (B.86) 



which can be substituted in the second gap equation (|B.85|) . yielding an equation for w% which 
does not contain u>2 anymore. The nominator of eq. (jB.86|) is different from zero, since filling in 
the numbers gives 

n2 ~ 1 ^ N 2 -lf 975 5984 \ ,„ . 

where we kept in mind that for a meaningful result, x ~ p-, should be positive. 

A numerical investigation of the gap equation (|B.85|) using ea. ()B.86|) revealed that there are 
no zeros. We conclude that there are no solutions with m 4 x~ 2a > A 4 x~ 2b . 

Next, let us find out if a possible solution with fh 4: x~ 2a = \ 4 x~ 2b or u)\ = might exist. 
We explicitly evaluate the gap equations (|B.77|) - (|B.78|) . where now Ti 2 = u\, 

N N 2 - 1 , . 3(N 2 - 1) u% uj 2 , 

W 32"7T Z 32"7T Z Wl — Ll>2 U>1 

^From the foregoing expressions, we infer that the limit lo\ —* ll>2 exists, giving rise to 

2(iV 2 - 1) 3(iV 2 - 1) 
2Cotj2 - 32tt2 + 32vr 2 - °' (R90) 

N 2 - 1 iV 2 - 1 3(iV 2 - 1) .„ N 

2 Ui ^- + -i = . (B.91 

iV 32tt 2 32tt 2 v 7 

This means that we have two equations to solve for the single quantity ui, which is present in 
I3i and &2 through the quantity x. It would be an extreme coincidence if these two different 
equations, which can be rewritten as 

- -3^' < R92 > 

111 - ~- (B.93) 
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Wi=W2 



1.2 

1 




Figure 5: The solution u\ = 0J2 as a function of 60 of ea. (|B.92|) . top curve, and ea, ()B,93j) . bottom 
curve, in units A^g = 1. Clearly, these two curves do no coincide. 



possess a common solution. That this is not the case can be inferred from the numerical solutions 
of both equations (|B.92|) and ()B.93|) . shown in Figure 5. 

As a final step, we should investigate if there is a solution in the region m 4 x _2a < X A x~ 2b . We 
can still define the coordinates uj\ and 102 by 



uJ 2 



rh 2 x a + iy —fh A x~ 2a + X 4 x~ 2b 



m 



2 x -a _ iy / _ff l 4 x -2a + \* x -2b 



(B.94) 
(B.95) 



In this case, u>i and UJ2 are complex conjugate. Henceforth, it would be more appropriate to use 
the modulus R and the argument <j>, <j> G] — vr, n], defined by 



Re* 
Re'* 



^2 , 



(B.96) 
(B.97) 



If the argument </> is so that \(f)\ > ^, then m 2 x a < 0. As a consequence, the estimate for (vl 2 ) 
will be positive. 

Most of the foregoing analysis can be repeated. The action (|B.81|) is rewritten in terms of 
R and 4> by 



(N 2 - 1) 

1 \ / 



2^ — tt-^UiR 2 + 2( U 2 R 2 cos 2 . 



N 

3R 2 (N 2 - 1 
32^2 



cos(2^>) ( In — 2 



R 5 
Ti 2 ~ 6 



>sin(2c 



The gap equations (|B.84|) - (|B.85|) reduce to 
TV 2 — 1 

-UtRe'* + Q U 2 R{e* + e~*) + 



N 



3(N 2 - l)Re* ( 1 
32^ 



+ k 



(B.98) 



(B.99) 



and its complex conjugate. With the parametrization (|B.96|) . we have fi 2 = R. 
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We must solve the following two real equations 10 for and R. 

- 2 ^\T^i cos <t> + 2 Co^2 cos + 3( " N 32 ~ 2 ~ (~ ~ ^ sin = 0, (B.100) 

iV 2 — 1 3(N 2 - 1) / sine* ,\ 
2^^isin0 + V 32 ^ 2 ; ( -^ + 0cos0j = 0. (B.101) 

We can divide these equations 11 by cos to obtain 

N 2 -1 3(N 2 - 1) ( 1 \ 

- 2— + 2C U 2 + V 32?r2 M ---0tan0j = 0, (B.102) 

iV 2 — 1 3(iV 2 -l)/tani \ 

2^^ 1 tan0 + ^ + = 0. (B.103) 

These equations can also be decoupled. The most efficient way to proceed is to eliminate R 
between these two equations to obtain an equation for 0, as the range in we must search for a 
solution is limited for this angle. The equation for finally becomes 

-909857V - 107712tt 2 6 + 12JV (4840 cot + 1734 In / -ii7(5O+ii^csc0sec0) \ _ 1573 ^ tan 



107712vr 2 

while the value of R is obtained from 



= 
(B.104) 



_ . R , 1950 + 4290 esc sec Ar .„ 

x = B In -p5- + b = „ —N . (B.105) 

pu A 2 — 11968^ 2 V ; 

MS 

We shall concentrate on the case N = 3. Depending on the value of the parameter bo, there is 
more than one solution possible. In Figure 6, we have plotted the expression HB.104j) for several 
values of the parameter b , namely b = 0.25,0, -0.25, -0.3, -0.33564..., -0.41594..., -0.5. It is 
possible to obtain those values of bo where the number of solutions change. If we consider the 
plots of Figure 6, it is apparent that for each 6 0) the corresponding curve possesses two extremal 
values. The number of solution exactly changes at those values of &o where the curve becomes 
tangent to the 0-axis. An explicit evaluation learns that his occurs at bo = —0.41595..., where 
= 2.26407... and at bo = —0.33564... where = 2.62545. It is important to know these num- 
bers to a high enough accuracy, to instruct the computer in which 0-interval it can search for a 
solution. If the initial values are not chosen in an appropriate way, the iterations will jump be- 
tween the different branches of solutions and there will be no convergence to any of them. There 
is a single solution if 6 > -0.33564... or b < -0.41595.... If -0.41595... < b < -0.33564..., 
there are three solutions, while for bo = —0.41595... and bo = —0.33564... there are two solutions. 
In Figure 7, we have displayed the solution for and R. To determine the solution which 
characterizes the vacuum, we should take that one which gives us the absolute minimum of the 
energy functional T, which was shown in Figure 3. 

As a final remark, we would like to notice that the same decomposition as in eq. (|4.7jl could 
also be useful for higher loop computations. The effective action T will remain symmetric under 
the exchange of uj\ and u>2 and equations like (|B.77|) - (jB.80|) shall remain valid. This should 
facilitate at least a bit the two-loop evaluation of the effective action and gap equations. Also, 
one does not need to evaluate any new anomalous dimension, since these are already known, 
either from previous calculations [31 1121 12*5] . or from exploiting relations like ea. ()B.41|) . 

10 The independence is hidden in LSi and U2 

11 We may assume coscj) 7^ 0, otherwise ecms. IB. 1001 - IB. 1011 would give <j> = 0, which is inconsistent with 
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gap eq. for phi 




phi 



Figure 6: The gap equation (|B. 104|) with N = 3 plotted in function of for the values b$ = 
0.25, 0, -0.25, -0.3, -0.33564..., -0.41594..., -0.5 (from bottom to top). 
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